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The multiplicity difference correlators between two well-separated bins in high-energy heavy-ion 
collisions are studied as a means to detect evidence of a first-order quark-hadron phase transition. 
Analytical expressions for the scaled factorial moments of multiplicity difference distribution are 
obtained for small bin size 5 with mean multiplicity in the bin s < 1.0 within Ginzburg-Landau 
description. The scaling behaviors between the moments are still valid, though they behave com- 
pletely different from the so-called intermittency patterns. A universal exponent 7 = 1.4066 is given 
" to describe the dynamical fluctuations in the phase transition in small 5 limit. 
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' The theoretical study of fluctuations in quark-hadron phase transition has been suggested for several years The 
^ > first motivation came from the quantitatively different results in Monte Carlo simulations on intermittency for pp 
^ ' collisions without phase transition from those theoretical predictions with the onset of phase transition^^l . Since then, 
multiplicity fluctuations are studied with phase transitions of second-order!^] and first-orderl"*"^] within Ginzburg- 
Landau model, and it is regarded as a possible means to reveal some features for the phase transition. Most of these 
works give the violation of the intermittency patterns but show remarkable scaling behaviors between Fq and 
and there seems to exist a universal exponent v jt jg suggested that the exponent v can be used as a useful 

diagnostic tool to detect the formation of QGP. In [4, 5] InFg are studied analytically for first-order phase transition 
and are expanded as power series of (J^/'^, and it is shown that the set of experimentally fitted coefficients for 5^^^ 
^vgj ] term can be used as a criterion for the onset and the order of the phase transition. 

■ It is known for a long time that the investigation of multiplicity fluctuations is very different in heavy-ion collisions, 
T— I \ though the power-law dependence of InF, on 5, Fq oc (5^'^«, has been found ubiquitous in hadronic and leptonic 

■ processes!^!. The main differences between heavy-ion physics and hadronic & leptonic ones on multiplicity fluctuations 
\ were noticed earlier in Ref. [8]. In Ref. [9, 10] an alternative way was proposed to study the fluctuations by means 
. of factorial moments of the multiplicity difference (FMMD) between two well-separated bins. This alternative is a 
' hybrid of the usual factorial correlators [-'^^1 and wavelets '"'^^l because Wjk in Haar wavelet analysis is just the difference 

of multiplicities in two nearest bins. When discussing the multiplicity difference, the two bins are not necessary the 
O • nearest ones. Instead, it will be assumed in present study that the two bins are well-separated. Let the two bins, each 
D ' of size 5"^ and separated by A, have multiplicities ni and 712, and define their multiplicity difference m — \ni — n2\. 

Let Qm be the distribution of multiplicity difference, which may be dependent on A, 5 and details of the process. 
^ ' Scaled FMMD are defined as 
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Moments defined above are similar to but not the same as the Bialas-Peschanski correlators [^^^ Fqiq2J because of (1) 
J-q are moments of the multiplicity difference between two bins and (2) J-q may depend on both A and 5 while Fq-^^q^ 
depends only on A. 

In Ref. [9], Tq are numerically studied within Ginzburg-Landau model. The scaling behaviors between Tq and 

T2, Tq (X T2'', are shown with Pq = [q — 1)^ and a universal exponent 7=1.099. In Ref. [10], J^q are analytically 
studied for very small 6 within the same model for second-order phase transition. The dynamical components of the 
moments are introduced and are shown to have similar scaling behaviors. The universal exponent 7=1.3424 given 
in [10] for dynamical fiuctuations is shown to be different from that in [9] but close to v given in the study of usual 
factorial moments I'^'^l. The closeness of 7=1.3424 to ~ 1.30 is reasonable because they both describe dynamical 
fiuctuations in phase transition. The little difference comes from the different x regions concerned: v corresponds to 
X region around — In a; ~ 1 but 7 to a; ^ 0, as discussed in [10]. The most important feature about this 7 is that it is 
completely determined by the general features of the model and independent of the parameters for the model. 

In this Letter, scaled FMMD Tq in first-order phase transition will be studied analytically for very small bin size 
d. For simplicity, the discussion is limited to two identical small bins under the condition that the mean multiplicity 
in each bin is less than or equal to 1.0. 

As a base and starting point, discuss the trivial and simplest case in which there is no correlations between the two 
bins and within each bin. Let the mean multiplicity in each bin is s. Because no dynamical reason is assumed, the 
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multiplicity distribution for each bin is a Poisson one 



Puis) = —r CXp(-s) . 

n] 

Prom this distribution, one can deduce the multiplicity difference distribution as 

Pm(s) = /m(2s)e-2^ {2-6m0) , 

where (2;) is the modified Bessel function of order m, 

(^/2)2'=+™ 



k=0 



k\{k + m)\ ' 



FMMD for pure statistical fluctuations are 



_^(stat) ^ ^ ^(^ _ 1) . . . _ ^ + l)p^(5) . 
m>q 

To complete the summation in last equation, one can introduce a generating functionl^l 

d''G{x,s) 



G{x, s) = 2e-2« x'"Imi2s) , Gq{x, s) = 
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With this function, y^^*'^*) can be rewritten as 



Direct algebra shows that 

G(a;, s) = 2e(^-2)' 
with a° = {-iys^y{ilf for i = 0, 1, • • •, and that 



/f-*)=G,(l,s) = Gg(s) . 

g 1 — exp(— xs) 
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Gg{s) = 2e- 



i=l j 
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where a| can be calculated by recurrence relation from a°, Oq = sOq"^, af = saf~^, af = saf~^ + afZi , (i > 2). 

In this Letter, we are interested only in small bin analysis, for which the mean multiplicity is less than 1. Then 
Gq{s) can be approximately written as 



Gg{s) = 2e-*s« 



1 + 
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(g + l)((7 + 2) (q + l)(g + 2)((7 + 3) (g + + 2)((/ + 3)((/ + 4) 



(9) 



and all other terms can be omitted. Direct estimation shows that the error caused by this approximation is less than 
1% for s < 1.0. 

For cases with dynamical fluctuations due to phase transition, the distribution of multiplicity difference is'^l 
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Qm{S,T)=Z-' / 2?</.P„(<5V I </> |^)e 



2\-F\, 



(10) 



where t is an indication of lifetime of the whole parton system, VcJ) = ird \ (/) l"^ , Z = J D^e ^^'^^ and the free energy 
F[(j)] = dz [a \ (f) +b \ (p +c \ (j) \^] for first-order phase transition. 
Substituting Qm{S,T) into Eq. (1), one gets 
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fg= dyGg(™a;2y)e-^'+"^'^+^^y dye-y"+ 
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with X = —b{5/c)'^/^ related with the bin width 5, u =\ a \ c/b"^. As discussed in Rcf. [4-6], b is negative for first-order 
phase transition. So x and u are both positive in present discussions. For the convenience of analytical calculations, 
define W 



/•oo 

Hg{u,v) = / y'' exp{-y^ + uv'^y + vy"^) 
Jo 
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which satisfies recurrence relations 

H2iu,v) = 1 + ^ {uv^Ho{u,v) + 2vHi{u,v)) , 
Hq+3{u, v) = ^ [{q + l)Hq{u, v) + uv'^Hq+i{u, v) + 2vHg+2{u, v)] . 

With Hq{u,v), fq can be expressed as 

(Tux^Hg+i 2{Tux^fHq+5 , 3{Tux^fHq+e 



Ho{u,x) 



Hq + r .ir^^:. - , :\:r : ::r . + 



(g + l)(g + 2) (g + l)(5 + 2)(g + 3) (g + + 2)(g + 3)(g + 4) 



(13) 



all Hq in the bracket in last equation should be read as Hq{—{T — l)w, x). 

The scaled FMMD obtained contain contributions from statistical fluctuations, contrary to the usual scaled 
factorial ones. One can see this clearly if one notices the fact that J-g are not equal to 1 for the case with pure 
statistical fluctuations. To eliminate the pure fluctuations from the moments, one can define the dynamical scaled 
FMMD as[i°l 

■^1 .j-(stat) ■ ^ ' 

J- q 

To make the definition sense, one should ensure that the mean multiplicity is the same for all the calculation of 
the moments concerned. In Ginzburg-Landau model, the mean multiplicity is s = tux"^ Hi{u,x)/ Hq{u,x). Then 

deviations of .Fg*^^"^ from one should indicate the existence of dynamical fluctuations. Different from the case for 

second-order phase transition, s oc in small x region for first-order phase transition, thus the pure statistical 
fiuctuations have no contribution to the slopes of InTg in the x region we arc now interested in. Thus, we will only 
discuss In J-g in present discussion. 

The dependences of InTg on — Inx from 2.0 to 4.0 are shown in Fig. 1 for t=10.0, u = 1.0 and 5.0. The — Ina; 
range is chosen from the requirement that s is much less than 1.0 for r = 10.0 and u = 5.0. Prom this figure, one can 
see clearly that In Fg depend strongly on parameters chosen, and In Tg decrease quickly with the increase of — In a; 
in the x range chosen. For sufficiently large — In a;, InTg satuate to values independent of u. In fact, the satuation 
values are independent of any parameter in the model. 

Scaling behaviors between Tg are shown in Fig. 2 for the same choices of parameters as in Fig. 1. One can see 

that the scaling behavior, Tg oc i is valid for either set of the parameters, and the slopes show weak dependence 
Because of the linearity of curves in Fig. 2 for Tq, Pq can be determined accurately. One can guess that for small x, 
Pg may be independent of all parameters in the model, and it is indeed the case. For very small x, one can get, using 

Hg{u,v)c.i[r{^) + vrii + i)] 

^(dyn) _ (g ~ l)ao,i/ao,o + a,,.i/ aq,Q - gai,i/cn,o 
' ao,i/ao,o + a2,i/a2,o - 2ai,i/ai,o 

with ttq^o = r(2^) and ttq^i = T{1 + |). f3q as function of ln(g — 1) is shown in Fig. 3. A perfect scaling behavior is 
shown 



Pq = {q-iy (16) 

with 7—1.4066. This exponent corresponds to the limit x ^ 0. In real experimental analysis, one fits curves for InJ^g 
and In at finite x, so that one should get an exponent with a little difference from the one given here. In fact, the 
experimentally obtained exponent should be less than 1.4066. 

It should be pointed out that the exponent can be obtained simply from the leading term in Gg by dropping off all 
non-leading terms because they are related to higher orders of x and have no contribution to 7 which is connected with 
properties of the moments in the limit x 0. The exponent given here is different from that in [10] for second-order 
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phase transition. The main reason is due to the introduction of | |^ term in the free-energy. One can see that 
1 (j) 1^ term is necessary for the study of first-order phase transition. In the case of second-order phae transition, the 
same term may also play a role. Only when —x ^ 1 can the —y^ term be neglected in Eq. (11), and one returns to 
the Ginzburg-Landau description of the phase transition studied in [3] . So the discussions in this Letter are different 
from those in former studies. In fact, if | </> |^ term plays an important role, 7 should be close to 1.4066 for both 
first-order and second-order phase transitions in small S analyses. With the increase of the importance of | </) |^ term, 
the exponent 7 for second-order phase transition can undergo an increase from 1.3424 in [10] to 1.4066 given here. 
But if there is no phase transition, the exponent 7 should be much less than 1.3424. Thus if an exponent 7 for the 
dynamical fluctuations is found near to 1.4, then the occurrence of quark-hadron phase transition can be pronounced. 

In summary, scaled FMMD are studied analytically within Ginzburg-Landau model in a kinetical region with mean 
multiplicity in single bin less than 1.0 for first-order quark-hadron phase transition. The dynamical fluctuations in 
FMMD are extracted, which give the same physical contents as the usual scaled factorial moments. Scaling behaviors 
between scaled FMMD are shown, and a truly universal exponent 7 = 1.4066 is given. 
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Figure Captions 

Fig. 1 Dependences of In^^ on — Ina; for r = 10.0, u = 1.0 and u = 5.0. From lower to upper are curves for 
g=2,3,4,5,6,7,8,9, respectively. 

Fig. 2 Scaling behaviors of InJ^q vs lnjF2 for the same choices of parameters as in Fig.l. From lower to upper are 
curves for 9=3,4,5,6,7,8,9, respectively. 

Fig. 3 Scaling behavior of In f3q vs ln{q — 1). 
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